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Abstract. A cartoon of the Effective Field Theory of many nucleon systems is drawn, 
concentrating on Compton scattering in the two nucleon system, and on nd scattering 
in the three body system. 



The purpose of this presentation is to give a concise introduction into the Effective 
Field Theory (EFT, for a review see e.g. [1]) of two and three nucleon systems 



as it emerged in the last three years. However, I can only give a "teaser" with 
a lot of words and figures and a few cheats in details, referring to the literature, 
esp. to the excellent proceedings of the INT-Caltech Workshops 1998 and 1999 [2]. I 
concentrate on work undertaken with J.-W. Chen, R.P. Springer and M.J. Savage in 
[3,4], P.F. Bedaque in [5,6], and F. Gabbiani in [6]. M. Birse's talk at this Workshop 
provides a more formal investigation of the EFT of the two nucleon system, and 
U.-G. Meifiner's alternative approach he presented here follows Weinberg's original 
suggestion [7] but needs to be studied further. 



Effective Field Theory methods are largely used in many branches of physics where 
a separation of scales exists. In low energy nuclear systems, the two well sepa- 
rated scales are, on one side, the low scales of the typical momentum of the process 
considered and the pion mass, and on the other side the higher scales associated 
with chiral symmetry and confinement. This separation of scales was explored with 
great success in the mesonic sector (Chiral Perturbation Theory [8]) and in the one 
baryon sector (Heavy Baryon Chiral Perturbation Theory [9]), producing a low en- 
ergy expansion of a variety of observables (see also the Chiral Perturbation Theory 
section of this Workshop). It provided for the first time a description of strongly in- 
teracting particles which is systematic, rigorous and model independent (meaning, 
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independent of assumptions about the non-perturbative QCD dynamics). 



Three main ingredients enter the construction of an EFT: The Lagrangean, the 
power counting and a regularisation scheme. First, the relevant degrees of freedom 
have to be identified. In his original suggestion how to extend EFT methods to 
systems containing two or more nucleons, Weinberg [7] noticed that below the A 
production scale, only nucleons and pions need to be retained as the infrared rele- 
vant degrees of freedom of low energy QCD. Because at these scales the momenta 
of the nucleons are small compared to their rest mass, the theory becomes non- 
relativistic at leading order in the velocity expansion, with relativistic corrections 
systematically included at higher orders. The most general chirally (and iso-spin) 
invariant Lagrangean consists hence of contact interactions between non-relativistic 
nucleons, and between nucleons and pions, with the first few terms of the form 
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where N = (^j is the nucleon doublet of two-component spinors and P % is the 
projector onto the iso-scalar- vector channel, P 1 ^ = -^(o^c 1 )^ 7 ^)*- ° ( r ) are the 
Pauli matrices acting in spin (iso-spin) space. The iso-vector-scalar part of the 
NN Lagrangean introduces more constants C% and interactions and has not been 
displayed for convenience. The field £ describes the pion, = v^E = e in /^, 
f n = 130 MeV. is the chiral covariant derivative = + V^, and the 
vector and axial currents are = + £^m0> A* = ^i&ntf ~ £^0- The 

interactions involving pions are severely restricted by chiral invariance. As such, 
the theory is an extension to the many nucleon system of Chiral Perturbation 
Theory and Heavy Baryon Chiral Perturbation Theory. Like in its cousins, all 
short distance physics - branes and strings, quarks and gluons, resonances like the 
A or a - is integrated out into the coefficients of the low energy Lagrangean. In 
principle, these constants could be derived by solving QCD or via models of the 
short distance physics like resonance saturation. The most common and practical 
way to determine those constants, though, is by fitting them to experiment. 

The EFT with pions integrated out (formally, g\ — in (1)) is valid below the 
pion cut and was recently pushed to very high orders in the two-nucleon sector [10] 
where accuracies of the order of 1% were obtained. It can be viewed as a sys- 
tematisation of Effective Range Theory with the inclusion of relativistic and short 
distance effects traditionally left out in that approach. 

Because the Lagrangean (1) consists of infinitely many terms only restricted by 
symmetry, an EFT may at first sight suffer from lack of predictive power. Indeed, as 
the second part of its formulation, predictive power is ensured only by establishing 



a power counting scheme, i.e. a way to determine at which order in a momentum 
expansion different contributions will appear, and keeping only and all the terms 
up to a given order. The dimensionless, small parameter on which the expansion 
is based is the typical momentum Q of the process in units of the scale A at which 
the theory is expected to break down, with estimates ranging from A w pa 300 
to 800 MeV [2] in the two body system for the theory with pions. The pion- 
less theory should be in disagreement with experiment starting at the pion cut, 
A n07r pa 140 MeV. Values for A and Q have to be determined from comparison 
to experiments and are a priori unknown. Assuming that all contributions are of 
natural size, i.e. ordered by powers of Q, the systematic power counting ensures 
that the sum of all terms left out when calculating to a certain order in Q is smaller 
than the last order retained, allowing for an error estimate of the final result. 

Even if calculations of nuclear properties were possible starting from the under- 
lying QCD Lagrangean, EFT simplifies the problem considerably by factorising it 
into a short distance part (subsumed into the coefficient of the Lagrangean) and 
a long distance part which contains the infrared-relevant physics and is dealt with 
by EFT methods. EFT provides an answer of finite accuracy because higher order 
corrections are systematically calculable and suppressed in powers of Q. Hence, the 
power counting allows for an error estimate of the final result, with the natural size 
of all neglected terms known to be of higher order. Relativistic effects, chiral dy- 
namics and external currents are included systematically, and extensions to include 
e.g. parity violating effects are straightforward. Gauged interactions and exchange 
currents are unambiguous. Results obtained with EFT are easily dissected for the 
relative importance of the various terms. Because only S'-matrix elements between 
on-shell states are observables, ambiguities nesting in "off-shell effects" are absent. 
On the other hand, because only symmetry considerations enter the construction 
of the Lagrangean, EFTs are less restrictive as no assumption about the underlying 
QCD dynamics is incorporated. 

In systems involving two or more nucleons, establishing such a power counting 
is complicated by the fact that unnaturally large scales have to be accommodated, 
so that some coefficients in the Lagrangean may not be of natural size and hence 
possibly jeopardise power counting: Given that the typical low energy scale in 
the problem should be the mass of the pion as the lightest particle emerging from 
QCD, fine tuning seems to be required to produce the large scattering lengths in 
the % and 3 Si channels (l/a ls ° = -8.3 MeV, l/a 3Sl = 36 MeV). Since there is a 
bound state in the 3 Si channel with a binding energy B = 2.225 MeV and hence a 
typical binding momentum 7 = \JMB ~ 46 MeV well below the scale A at which 
the theory should break down, it is also clear that at least some processes have to 
be treated non-perturbatively in order to accommodate the deuteron. Most likely, 
these small scales do not arise from the fact that the real world is close to the chiral 
limit: In the singlet channel, for instance, the one pion exchange potential vanishes 
in the chiral limit and thus cannot be the cause of the fine tuning. The fine tuning 
then must be a result of short distance physics. 

A way to incorporate this fine tuning into the power counting was suggested by 



Kaplan, Savage and Wise [11]: At very low momenta, contact interactions with 
several derivatives - like p 2 C*2 and the pion-nucleon interactions - should become 
unimportant, and we are left only with the contact interactions proportional to Cq. 
The leading order contribution to nucleons scattering in an S wave comes hence 
from four nucleon contact interactions and is summed geometrically as in Fig. 1 to 
all orders to produce the shallow real bound state, i.e. the deuteron. 

= ^X + M + )oo( + - = T f^ 

FIGURE 1. Rc-summation of the contact interactions into the deuteron propagator. 



How to justify this? Any diagram can be estimated by scaling momenta by a 
factor of Q and non-relativistic kinetic energies by a factor of Q 2 /M. The remaining 
integral includes no dimensions and is taken to be of the order Q° and of natural 
size. This scaling implies the rule that nucleon propagators contribute one power 
of M/Q 2 and each loop a power of Q 5 /M. Assuming that 

c ° ~ Wq ' ° 2 ~ mXq 2 ' (2) 

the diagrams contributing at leading order to the deuteron propagator are indeed 
an infinite number as shown in Fig. 1, each one of the order 1/(MQ). The regu- 
lator dependent, linear divergence in each of the bubble diagrams does not show 
in dimensional regularisation as a pole in 4 dimensions, but it does appear as a 
pole in 3 dimensions which we subtract following the Power Divergence Subtrac- 
tion scheme [11]. Dimensional regularisation is chosen to explicitly preserve the 
systematic power counting as well as all symmetries (esp. chiral invariance) at each 
order in every step of the calculation. At leading (LO), next-to-leading order (NLO) 
and often even NNLO in the two nucleon system, it also allows for simple, closed 
answers whose analytic structure is readily asserted. The deuteron propagator 

in (3) 

has the correct pole position and cut structure when one chooses 

Indeed, when choosing // ~ Q, the leading order contact interaction scales as de- 
manded in (2) and - as expected for a physical observable - the NN scattering 
amplitude becomes independent of /i, the renormalisation scale or cut-off chosen. 
The same can be shown for the higher order coefficients, so that the scheme is 



self-consistent. Power Divergence Subtraction moves hence a somewhat arbitrary 
amount of the short distance contributions from loops to counterterms and makes 
precise cancellations manifest which arise from fine tuning. Notice that the re- 
summed deuteron propagator has the same order 1/(MQ) as each diagram in Fig. 1. 

One surprising result arises from this analysis because chiral symmetry implies 
a derivative coupling of the pion to the nucleon at leading order. The contribution 
from one pion exchange includes a factor of Q~ 2 from the pion propagator and a 
factor of Q 2 coming from the pion-nucleon vertices, so that for momenta of the 
order of the pion mass, the instantaneous one pion exchange scales as Q° and is 
smaller than the contact piece Co which according to (2) scales as Q^ 1 . Iterated and 
radiative pion exchanges are suppressed even further. Pion exchange and higher 
derivative contact terms appear hence only as perturbations at higher orders. In 
contradistinction to iterative potential model approaches, each higher order con- 
tribution is inserted only once. In this scheme, the only non-perturbative physics 
responsible for nuclear binding is extremely simple, and the more complicated pion 
contributions are at each order given by a finite number of diagrams. For example, 
the NLO contributions to the deuteron are the one instantaneous pion exchange 
and the four nucleon interaction with two derivatives, Fig. 2. The constants are 
determined e.g. by demanding the correct deuteron pole position and residue [12]. 
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FIGURE 2. The NLO corrections to the deuteron. 



In the two body sector, the theory thus emerging has been put to extensive tests at 
NLO and NNLO, giving for the first time analytic answers to many deuteron prop- 
erties, see e.g. [2]. Although in general process dependent, the expansion parameter 
is found to be of the order of | in most applications, so that NLO calculations can 
be expected to be accurate to about 10%, and NNLO calculations to about 4%. In 
all cases, experimental agreement is within the estimated theoretical uncertainties, 
and in some cases, previously unknown counterterms could be determined. 

The elastic deuteron Compton scattering diagrams to NLO are partially obtained 
by gauging the Lagrangean (1), i.e. by replacing ordinary derivatives by covariant 
ones: At LO, a seagull-graph and one graph in which the incident and outgoing 
photon couple to the same nucleon are found. At NLO, the photons are attached in 
all possible ways to the corrections in Fig. 2, including to the pion, the NNir vertex 
and the C 2 vertex. The Fermi interaction & ■ B probes the 1 S intermediate NN 
state and enters at NLO, too. Finally, the iso-scalar electric nuclear polarisability 
was shown to come from relativistic ( "radiative" ) pions in Chiral Perturbation The- 
ory [13], Q!e : n = 4Snpm ' anc ^ * s NLO. The cross section fits finally on less than one 
page with functions not more complicated than Logarithms and Arcustangentes [4] 



and contains no free parameters. Comparison with the Urbana experiment [14] in 
Fig. 3 shows good agreement, with the pion graphs that dominate the electric po- 
larisability of the nucleon necessary to improve it. The deuteron scalar and tensor 
electric and magnetic polarisabilities are also easily extracted [3]. 



30 

da 
aft 

20 

(nb/sr) 

10 





E 1 = 49 MeV 








50 



100 



150 

( degrees ) 



E 1 = 69 MeV 




50 



100 150 

8 ( degrees ) 



FIGURE 3. The differential cross section for elastic 7-deuteron Compton scattering at incident 
photon energies of E~ = 49 MeV and 69 MeV in an EFT with explicit pions [4], no free parameters. 
Dashed: LO; dotted: NLO without the graphs that contribute to the nucleon polarisability; solid 
curve: complete NLO result. 



In the three body sector, even the leading order calculation is too complex for a fully 
analytical solution. Still, the equations that need to be solved are computationally 
trivial and can furthermore be improved systematically by higher order corrections 
that involve only (partly analytical, partly numerical) integrations, as opposed to 
many-dimensional integral equations arising in other approaches. The nd system 
provides a laboratory in which many complications of the other channels are not 
encountered: The absence of Coulomb interactions ensures that only properties 
of the strong interactions are probed. In the quartet channel, the Pauli principle 
forbids three body forces [15] in the first few orders. Because the calculation is 
parameter- free, it allows one to determine the range of validity of the KSW scheme 
without a detailed analysis of the fitting procedure. Although e.g. the quartet 
scattering length is large, no extra fine tuning except the one for the deuteron 
is required. In the S wave, spin-doublet (triton) channel, the situation is more 
complicated. An unusual renormalisation of the three-body force makes it large 
and as important as the leading two-body forces [16]. More work is needed there. 

A comparative study between the theory with explicit pions and the one with 
pions integrated out was performed in [5] for the spin quartet S wave. As seen above, 
the two theories are identical at LO: All graphs involving only Co interactions are 
of the same order and form a double series which is not geometrical and cannot 
be summed analytically. One is hence left with the task of summing all "pinball" 
diagrams (first line of Fig. 4). Summing all "bubble-chain" sub-graphs into the 
deuteron propagator, one can however obtain the solution numerically from the 
integral equation pictorially shown in the lower line of Fig 4. A code runs within 
seconds on a personal computer. 

The power counting shows that at NLO, we have additional contributions from: 




FIGURE 4. The double infinite series of LO "pinball" diagrams, some of which are shown in 
the first line, is equivalent to the solution of the Faddeev equation shown in the second line. 



p 2 C*2 insertions and pion exchange corrections to the deuteron propagator depicted 
in the first line of Fig. 5; pionic vertex corrections to Co (second line); and the 
pion diagram of the last line which corrects the three particle intermediate state. 
Here, we used a re-formulation of the Lagrangean (1) in order not to have poorly 
convergent diagrams containing C 2 like the second one in the first line of Fig. 4, 
see [5] for details. The calculation without explicit pions was carried out to NNLO. 
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FIGURE 5. The NLO contributions to nd scattering in the quartet channel. First line: Cor- 
rections to the deuteron propagator; second line: pionic corrections to the Co vertex ; third line: 
pionic corrections to three particle breakup in the intermediate state. Permuted graphs left out. 

All calculations demonstrate convergence. The scattering length is a( 4 «S| , LO) = 
(5.1 ± 1.5) fm, and at NLO with (without) perturbative pions a( 4 S , 3 , NLO, ir) = 
(6.8 ± 0.7) fm (a( 4 5|,NLO,novr) = (6.7 ± 0.7) fm). At NNLo' [15] report 
a( 4 Sp NNLO, no7r) = (6.33 ± 0.1) fm, and the experimental value is a( 4 S|, exp) = 
(6.35 ± 0.02) fm [17]. Comparing the NLO correction to the LO scattering length 
provides one with the familiar error estimate at NLO: (|) 2 ~ 10%. The NLO 



calculations with and without pions lie within each other's error bar. The NNLO 
calculation is inside the error ascertained to the NLO calculation and carries itself 
an error of about (|) 3 ~ 4%. NLO and LO contributions become comparable for 
momenta of more than 200 MeV. In the imaginary part shown in Fig. 6, the same 
pattern emerges with a slightly more pronounced difference between the pion-less 
and pions-full theory. Because results obtained with EFT are easily dissected for 
the relative importance of the various terms, one concludes that pionic corrections 
to nd scattering in the quartet S wave channel - although formally NLO - are 
indeed much weaker: The calculation with perturbative pions and with pions inte- 
grated out do not differ significantly over a wide range of momenta. The difference 
should appear for momenta of the order of and higher because of non-analytical 
contributions of the pion cut. However, it is very moderate for momenta of up to 
300 MeV in the centre-of-mass frame (E cm « 70 MeV), see Fig. 6. This and the 
lack of data makes it difficult to assess whether the KSW power counting scheme 
to include pions as perturbative increases the range of validity over the pion-less 
theory, but effects from the pion cut are seemingly weak. 
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FIGURE 6. Real and imaginary parts in the quartet S wave phase shift of nd scattering versus 
the centre-of-mass momentum [5]. Dashed: LO; solid (dot-dashed) line: NLO with perturbative 
pions (pions integrated out); dotted: NNLO without pions [15,6]. Realistic potential models: 
squares from [18], crosses from [19], triangles from [20]. Stars: TUNL pd phase shift analysis [18]. 

Finally, the real and imaginary parts of the higher partial waves / = 1, . . . , 4 in the 
spin quartet and doublet channel were presented in [6] in a papameter-free calcu- 
lation. Figure 7 shows two examples. Comparison of the LO with the NLO and 
NNLO result demonstrates convergence of the EFT, with the expansion parameter 
again about |. It is interesting that the NLO correction is for high enough energies 
sometimes sizeable, while the NNLO correction is in general very small. Within the 
range of validity of this pion-less theory, convergence is good, and the results agree 
with potential model calculations (as available) within the theoretical uncertainty. 
That makes one optimistic about carrying out higher order calculations of prob- 
lematic spin observables like the A y problem where the EFT approach will differ 
from potential model calculations due to the inclusion of three-body forces. 
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FIGURE 7. Real and imaginary parts of the quartet P (top) and doublet F (bottom) wave phase 
shift of nd scattering versus the centre-of-mass energy in the EFT without pions [6]. Dashed: LO; 
dot-dashed: NLO; solid line: NNLO. The experiments by Huttel et al. [21] and Schmelzbach et 
al. [22] are denoted by open squares and diamonds, respectively. The calculations of Kievsky et 
al. (crosses) are from Refs. [19] below breakup (E cm = B) and [23] above breakup. 
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